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We have performed fluorescence x-ray absorption fine-structure (XAFS) measurements from 
20-200 K on a 5000 A, c-axis film of YBa 2 Cu 3 7 _i (YBCO) on MgO (T e =89 K) using photons 
polarized perpendicular to the film. The quality of the data is high out to 16 A -1 . The data from 
3-15.5 A" 1 were transformed into r-space and fit to a sum of theoretical standards out to 4.0 A. 
These data are compared to YBCO data from a single crystal and from a film on LaAlC>3 with the 
same T c . The main difference between the single crystal and the film data is that while the single 
crystal data is well described by a two-site axial oxygen [0(4)] distribution, we see no evidence for 
such a distribution in either thin film sample. We place an upper limit on the size of an axial oxygen 
splitting for the film on MgO at A r <0.09 A. Therefore, the magnitude of the splitting is not directly 
related to T c . Fits to the temperature dependent data from the YBCO film on MgO indicate that 
all bonds show a smooth change of their broadening factor a, except the Cu-0(4) bonds, which show 
an increase in a in the vicinity of T c , followed by a decrease of the same magnitude. Such a feature 
does not occur in diffraction measurements. Since XAFS measurements of a include any correlation 
between the atoms in a given bond, we conclude that the 0(4) position becomes less correlated with 
the Cu positions near T c . Correlation measurements of these and several further near-neighbors are 
also reported. 
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I. INTRODUCTION 

Since the discovery of the high T c superconductors!] 
there have been several measurements of possible struc- 
tural changes near the superconducting transition. Some 
measurements were far enough above the noise that the 
result was immediately convincing to the scientific com- 
munity, such as Sharma et aFs ion channeling mea- 
surements of a sudden correlation of the alignment in 
the a6-plane_pf the Cu(l), 0(4) and Cu(2) .atoms in 
ErBa 2 Cu 3 7 cl and in YBaaCuaOr-a (YEGG).! (Fig. 
shows structure and site definitions for YBCO.) Some 
have relied on relatively small changes in a relatively large 
signal, such as Horn et aPta diffraction measurement of 
an anomaly in the orthorhombicity [(b — a)/(b + a)] of 
YBCO, or Mustre de Leon et aFs x-ray-absorption fine 
structure (XAFS) measurement of a split 0(4) position 
distribution (Sr ~0.13 A) in YBCO thatpde creases some- 
what or disappears (Sr <0.11 A) near T c J§eI Both of these 
measurements were initially treated with a fair amount of 
skepticism, because obtaining a reliable signal for mea- 
suring such small changes is difficult. Horn et-.aPs re- 
sult has since been confirmed both by-diffractionEI and by 
capacitive dilatometry measurementsu and consequently 
must be considered real, although its relation to T c is still 
unclear. On the other hand, Mustre de Leon et aFs mea= 
surement has only been partially confirmed by XAFS0 
(see below), and no other experiment has been able to 
verify the result in YBCO. 

The XAFS measurement of a split 0(4) site distribu- 
tion has been taken as a possible explanation for anomiu 
lous features in vibrational spectra involving the 0(4).t£l 
This distribution may also fits into current polaron mod- 
els for electron-phonon-coupling (for instance, see Refs. 
[lfL|l4]), In spite of these important implications, the pos- 
sibly dynamic nature of the 0(4) distribution has never 
been verified. In fact, it has been suggested that the split 
distribution is actually due to ordemd-P(l) vacancies ,Ej 
as in the Ortho-II phase of YBCO.BU Part of the am- 
biguity in the interpretation lies in the difficulty of the 
measurement itself. The most serious problems with-thc 
XAFS measurements are pointed out by Stern et aln In 
their attempt to confirm the Mustre de Leon et al mea- 
surements, they tried several samples with various oxy- 
gen concentrations. The only sample that was clearly 
consistent with a split distribution of the 0(4) position 
had the lowest T c of the samples measured (89 K). Fits 
to a split distribution for the other samples were less 
convincing; both the single-site fit and the split-site fit 
were unsatisfactory. The sample with the highest T c (92 
K) did not show any temperature dependence near T c , at 
least as far as the fits were concerned. Since the measure- 
ment relies on such a small change in the overall signal 
and the results were severely sample dependent, the au- 
thors concluded that an independent method is necessary 
to confirm the split-site character of the 0(4) distribu- 
tion. They also concluded that the relation between the 



split distribution and T c must be weak and unimportant 
to superconductivity since the highest T c sample showed 
the least change with temperature. This conclusion is 
weak, however, because the fits to this T c =92 K sample 
were of poor quality. 

Unfortunately, no independent experiments have been 
able to confirm the result in YBCO, although x-ray 
diffraction analysis has seen a.-(jiuch larger effect in sin- 
gle crystals of YBa 2 Cu2.7*07,li2 pulsed neutron diffrac- 
tion (using a PDF analysistj has seen correlated displace- 
ments of the Cu and axial oxygen in Tl-2212, andrXAFS 
has seen axial oxygen changes near T c in T1-1234.E3 

As Stern et al have pointed out, until an independent 
experiment confirms or refutes the possibly split nature 
of the 0(4) position, the split nature of the 0(4) site 
will be questionable. However, with different samples, 
analysis techniques, etc. XAFS can still help clarify the 
issue. For instance, well oriented samples are important 
for these XAFS measurements because the polarization 
of the incoming synchrotron light can be utilized to re- 
move the contribution to the XAFS of the in-plane oxy- 
gens. Magnetically-oriented powders were used in both 
the Mustre de Leon et al and the Stern et al studies. 
Each of these absorption measurements were made in 
the transmission mode. In the study presented here, we 
use a 5000 A film on MgO, a 4000 A film on LaA10 3 , 
and a 17 ^m (2% Ni) single crystal of YBCO, all with 
a T c of 89 K. Each of these samples are highly oriented 
with the c-axis perpendicular to the sample surface. All 
measurements were made in the fluorescence mode. Our 
analysis focuses on the film on MgO and uses the other 
two samples for comparison. 

In contrast to the results of Stern et al and Mustre 
de Leon et al, the film data can be well described by a 
single, harmonic 0(4) distribution that shows only very 
subtle changes near T c . This data is very reproducible 
over a wide range of temperatures, and thus provides a 
null result to compare with any other results. On the 
other hand, the single crystal data shows a split 0(4) 
distribution below T c consistent with previous measure- 
ments on oriented powders. By comparing a case where 
the splitting is not present and one where it is, we show 
that previous XAFS measurements of a split 0(4) distri- 
bution are not relying on an unusually small signal. 

The measurements presented here also provide details 
of the harmonic broadening of the pair distribution of 
the atom pairs in the film on MgO. These measure- 
ments show interesting temperature dependences which 
are interpreted as changes in the correlations between 
the atomic positions in a given atom pair. Fluctuations 
in the broadening parameters of the Cu(2)-0(4) pair, and 
to a lesser degree in the Cu(l)-0(4) pair in the vicinity 
of T c are also reported. 

We begin our report by describing the experimental 
details in Sec. [n]. The data analysis procedures are given 
in Sec. III. The isolation of the 0(4) contribution to the 
XAFS is presented in Sec. 



IV 



Results of detailed fits 
to the film data are given in Sec. [y| Consequences of 
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observed features in the data are discussed in Sec. VI 
and the conclusions of this work are summarized in Sec. 



VII. 



II. EXPERIMENTAL DETAILS 

The thin film samples of YBCO were prepared by B. 
Lairson, who used Molecular Beam Epitaxy (MBE) to 
deposit the superconductor onto a MgO and a LaAlC>3 
substrate. The film on MgO is estimated to be 5000 A 
thick, and the film on LaAlC>3 is approximately 4000 A 
thick. X-ray diffraction measurements indicate less than 
1% impurity phases. The diffraction measurements also 
indicate that the films are oriented with their c-axes per- 
pendicular to the plane of the film with less than 1% cl- 
ot 6-axis impurities for the film on MgO, and <3% for the 
film on LaAlOs. TEM pictures confirmed that films cre- 
ated in the same way are c-axis oriented, with trace im- 
purities and that the crystals are twinned. The samples 
have a c-axis lattice parameter of 11.64 A at room tem- 
perature. The transition temperatures were measured 
with the standard four-probe technique and found to be 
89 K. 

We also obtained a fairly large single crystal from W. 
Hardy which contains about 2% Ni. This crystal is also 
c-axis oriented. The thickness is estimated to be ~ 17 
/j,m by looking at the Cu K-edge step in the x-ray trans- 
mission spectrum. 

X-ray fluorescence measurements were made on beam- 
line 7-3 at the Stanford Synchrotron Radiation Labora- 
tory (SSRL). The samples were placed in an Oxford he- 
lium cryostat and oriented such that the incident x-ray 
beam was striking the surface of the samples between 10- 
12 degrees. Temperature was regulated via two sensors, 
one in the helium chamber and one about two inches 
above the sample. Because of the geometry and later 
calibration measurements, we estimate that the absolute 
temperature of the sample is probably 3-5 K higher than 
the sample sensor was indicating (the "nominal" temper- 
ature) . Relative changes in the temperature of these data 
should be much better and is probably around 1 K. All 
reported temperature measurements of the XAFS data 
are the nominal temperature. 

The fluorescence x-rays were detected using a 13- 
element Ge detector. This detector allows for energy res- 
olution that isolates the Cu K a peak and eliminates any 
need for background fluorescences to be removed from 
these data. Count rates in each channel were kept well 
below saturation levels (< 20% of saturation), but were 
in any case corrected for the dead-time-constant (~ 5- 
10 fj&) introduced by the energy-resolving peak-shaping 
filters. Besides the ability to resolve the Cu K a fluores- 
cence, the main utility in using the Ge detector is to aid 
in the removal of Bragg peaks in the absorption spec- 
trum. Because the Bragg spikes occur at the incident 
beam energy, the energy resolving power of the Ge detec- 



tors can partially discriminate against them. However, a 
large change in the fluorescence flux can bring the detec- 
tors much closer to saturation, and thus the dead-time 
correction becomes more important. If a Bragg spike is 
particularly strong, the dead-time correction cannot ac- 
count for the additional flux accurately enough, and the 
affected data must be removed. Such removal is accom- 
plished by replacing the affected data by the normalized 
mean of the other channels. In this way we can usu- 
ally generate a full absorption spectrum over the energies 
of interest with almost no experimentally induced spikes 
or dips. This method provides a significant advantage 
over more traditional single channel ionization counters 
when one realizes that for thin films on well-oriented sub- 
strates, or for single crystals, Bragg spikes can affect data 
over tens of eV. 

Unfortunately features which we call "glitches" can 
still arise from nonuniformities in the sample coupling 
to nonuniformities in the incident beam.Ej'E3 Since these 
occur in all the channels of the Ge detector, we must 
remove such features by fitting a low (2nd or 3rd) order 
polynomial through neighboring (unaffected) data points 
and replacing the data with the fit. Fortunately, glitches 
often only affect between one and four successive data 
points, and so this procedure does not significantly alter 
the XAFS. 

An example of both the film data which has been pro- 
cessed as described above is displayed in Fig. ||. The 
single crystal data is of similar quality. 



III. DATA ANALYSIS 

Our data analysis procedures have been detailed else- 
where for standard transmission experiments. In this 
work, we will describe briefly the general procedure and 
only give details where significant differences in the pro- 
cedure arise over Refs. 



A. Basic XAFS 

The main features in an x-ray absorption spectrum are 
due to single electron excitations, namely, the photoelec- 
tric effect. These excitations appear as sudden jumps in 
the absorption coefficient n(E). In a solid, the data just 
above these jumps exhibit oscillatory structure which has 
been named X-ray Absorption Fine Structure (XAFS). 
The XAFS are isolated by defining a function 



X(E) 



ME) 



(1) 



which is used just beyond the edge, normally > 25 eV. 
fio (E) is a normalization function which is the absorption 
due to all processes that do not include the photoelectron 
backscattering effect described presently. 
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The oscillations in /i(-B) are primarily due to an inter- 
ference effect between the outgoing photoelectron's wave 
function and the part of the photoelectron's wave func- 
tion which has scattered off nearby neighbors and re- 
turned to the absorbing atom. As the photoelectron's 
wave vector is increased, the interference is modulated 
with a frequency given by twice the distance to the 
backscattering atom. There are several other factors 
which must be taken into account. The version of the 
"XAFS equation" that we use in our data analysis is 

/>oo 

x{k)=ImJ2Ai 9l {r)e^ r+25 ^ +5 ^dr (2) 

where the amplitude factor Aj is given by 

A> = N^F^k), (3) 

Ni is the number of equivalent atoms in shell i, S 2 is 
the "amplitude reduction factor" that accounts primar- 
ily for many-body effects that reduce the XAFS oscilla- 
tion amplitude such as shake-up or shake-off, Fi(k) is the 
backscattering amplitude of the photoelectron off neigh- 
bors i including a reduction sue to the mean-free path 
of the photoelectron, gi(r) is the pair distribution func- 
tion for the absorbing and backscattering atoms, and the 
phase shifts of the photoelectron from the central and the 
backscattering atom are give by S c (k) and Si(k), respec- 
tively. The distribution function gi(r) is usually taken to 
be harmonic 

-(■—Hi) 2 

9i (r) = e~^r, (4) 

where Ri is the average distance between the central and 
the backscattering atoms i, and the broadening factor <ji 
(otherwise known as the Debye- Waller factor) combines 
the effects of thermal and static disorder. If a shell i 
contains atoms that are inequivalent, then a% may also 
be used as a measure of any distortion in the shell, even 
though such distortions are rarely harmonic. 

In studies of well-ordered materials (such as YBCO) 
we expect there to be very little difference in the bond 
lengths measured by XAFS compared to diffraction. In 
fact, with the exception of the 0(4) site, there are 
many XAFS studies which demonstrate the strong agree- 
ment between local structure in YBCO r^and, most of 
its relatives) and the average structure.Ea c§ On the 
other hand, comparative measurements of the correlated 
Debye- Waller factors given by XAFS and the uncorre- 
cted ones given by diffraction have been virtually ignored 
by most experimenters. 

Part of the reason for not making direct compar- 
isons between XAFS and diffraction broadening fac- 
tors is that getting the absolute broadening factor is 
not always possible. Experimental measurements of 
amplitude functions have their own (usually unknown) 
broadening, and until recently, theoretical calculations 
were not of high enough quality to give reliable results. 



FEFF6EM written by Rehr, Zabinsky and co-workers, 
has been shown to calculate accurate backscattering am- 
plitudes and phase shifts, and even gives reasonable 
values for Debye- Waller parameters far simple materi- 
als with known Debye temperatures. EStil Such calcula- 
tions provide a backscattering amplitude function with a 
known (zero) width for absolute comparison to real atom 
pairs. Although only experience will allow us to deter- 
mine how accurate absolute measurements of correlated 
broadening factors really are, our limited experience in- 
dicates that they are probably within 10% in most cases. 

Another reason direct comparisons are rarely made is 
simply that diffraction and XAFS width parameters mea- 
sure different quantities. Diffraction relies on coherent 
diffracting of x-rays off many crystal lattice planes to 
generate a diffraction spot, and therefore its measure- 
ments of atom positions and broadening factors gives an 
average over many unit cells. XAFS depends only on the 
instantaneous (within ^10~ 15 sec) position of neighbor- 
ing atoms with respect to the central (absorbing) atom, 
and therefore the broadening factors are a measure of 
variations in an atom-pair's distance. For example, if 
two neighboring atoms are vibrating with the same (spa- 
tial) amplitude W but are somehow rigidly connected, 
the bond length would never change, and XAFS would 
measure a zero width. Diffraction would measure broad- 
ening factors for the two atoms to be W. Mathematically, 
if 5rx is the instantaneous deviation of atom A's posi- 
tion from its mean position, the average deviation in the 
distance between atoms A and B is given by 

(T AB = \/< {Sta - Sr B ) 2 > 

= yj< Sr\> + < Sr 2 B > -2 < Sr A Sr B >■ (5) 

If we define a a = \J < Sr\ > (which equals y/U in con- 
ventional diffraction notation), then the last term can 
vary between — 1oa&b and +2&A &b, depending on the 
degree to which the motions of the two atoms are corre- 
lated. One can define a "correlation" parameter <j> from 

a AB = cr A + a B - 2vaCTb4>- (6) 

Since a\ is given by diffraction experiments, we can pro- 
vide measurements of correlation factors. <f> will vary 
between +1 for atoms vibrating in phase with each other 
(like in an acoustical phonon), for atom-pairs that are 
separated by a few unit cells at temperatures well above 
the Debye temperature, and —1 for pairs that are moving 
out-of-phase with each other (like in an optical phonon) . 
Of course since the XAFS measurements give average 
atom-pair distances over the crystal, static disorder can 
produce the same results. 

B. Data Reduction 

An initial estimate of x(&) is obtained by approximat- 
ing jJ-o(E) in Eq. |lj by a cubic spline though n{E), k by 
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^/2m(E E o)_ ^ ^ ^ e energy at the half-height of 
the edge. The data is fit in r-space as described in Ref. 
|23| . The fit is then subtracted from the data to obtain a 
residue which can be utilized to make a better estimate 
for fJ.o(E), as described in Ref. |2jj. 

The data have been corrected for the additional "self- 
absorption" of the fluorescing photon in the sample using 
a treatment which goes slightly beyond the usual correc- 
tion (for instance, see Troger et ao) by including the 
finite thickness of the sample and the effect of the XAFS 
oscillations in the correction term. The correction for the 
5000 A film is small (~ 8% of Ho(E)), but for the 17 fim 
single crystal the correction is more significant (~ 60% 
of ME))- 

An example of k 3 x(k) for both the film and the single 
crystal which has been reduced using the above proce- 
dures is displayed in Fig. ||. The Fourier transform (FT) 
of kx(k) for these data are shown in Fig. 0. 



IV. ISOLATING THE XAFS SIGNAL DUE TO 
THE 0(4) CONTRIBUTION 

Since the information above 2.5 A is complicated, we 
had to perform fits to the data to discern any informa- 
tion about the Y, Ba and Cu environment. These fits 
are described in Sec. [v| However, since the Cu(l)-0(4) 
and Cu(2)-0(4) peaks are relatively well separated, we 
can isolate them and back-transform the data to fc-space. 
Mustre de Leon et al based their argument that the 0(4) 
position was split by observing anharmonic behavior in 
the Fourier filtered fc-space data in the form of a beat 
near 12 A -1 . If a two-site distribution exists, we expect 
the contributions of the two atom-pairs to be split into 
four distinct distances: two corresponding to the mean 
Cu(l)-0(4) distance ± the splitting (i?c«(i)-o(4) =t A r ) 
and two corresponding to the mean Cu(2)-0(4) distance 
± the splitting (i?c«(2)-o(4) ± A r ). Since the XAFS 
signal for each of these distances is oscillatory, a small 
splitting will generate a beat in the XAFS signal at a 
A;- vector given by 2k A,. = it. In order to make our data 
directly comparable to their's, we have also isolated the 
0(4) contribution in this way. First, the FT of k 3 x(k) 
was obtained as in Fig. [|. We chose fc 3 -weighting for this 
transform because the 0(4) peak is much better resolved 
than in a fc-weighted transform. The Fourier transforms 
of the k 3 x(k) data were then back-transformed from 1.3 
to 2.2 A . The 0(4) contribution from this procedure is 
displayed for the single crystal in Fig. [s] and for the film 
on MgO in Fig. ^| The sinusoidal behavior in both these 
figures is the sum of the Cu(l)-0(4) and the Cu(2)-0(4) 
components. 

The 0(4) contribution in the single crystal displayed 
in Fig. |5j is in rough agreement with the original work 
of Mustre de LeonO; although the data in the 11-13 A -1 
range varies from temperature to temperature, the T=83 
K and the T=86 K data show clearly anharmonic behav- 



ior near 12 A -1 , indicating a splitting of about A r ~ 
0.13 A. In fact, except for the T=100 K data, none of 
the single crystal data can be fit wi th a single, harmonic 
distribution for the 0(4) site (Sec. VB ). 

The data from the film on MgO in Fig. ||, on the 
other hand, displays harmonic behavior at all tempera- 
tures and shows very little change from temperature to 
temperature. This harmonic behavior persists out to 16 
A -1 . Since a beat may still occur at higher wave vectors, 
these Fourier-filtered A-space data set an immediate up- 
per limit of 0.10 A. The character of the film data is 
distinctly different than the single crystal data, showing 
little if any splitting. This result also applies to the film 
on LaA103 and therefore appears to be a generic result 
of YBCO films. 



V. FITS 

Although the raw data has provided much useful in- 
formation about the nature of the 0(4) site distribution, 
we can still learn more details both about the 0(4) dis- 
tribution and about the further neighbors by performing 
fits of the spectra to the sum of theoretical standards. 
Such fits provide more quantitative details about the lo- 
cal environment around the copper atoms and allow us to 
deconvolve the contributions of overlapping peaks due to 
the Cu(2)-Y, Cu(2)-Cu(2) (between the planes), Cu(2)- 
Ba, Cu(l)-Ba and the Cu(l)-Cu(2) pairs. Furthermore, 
by fitting the data to theoretical standards calculated by 
FEFF6 we can obtain absolute estimates of the broaden- 
ing factors a and compare them to results from diffrac- 
tion studies, thus allowing for measurements of the cor- 
relation parameter tfi (Eq. ||). 

We only report fit results as a function of tempera- 
ture for the YBCO film on MgO data beyond the near- 
neighbor Cu(l)-0(4) and the Cu(2)-0(4) atom pairs. 



A. Fitting procedures 

Although weighting the XAFS x(k) by a factor of fc 3 
enhances the high k data where the beat in the 0(4) part 
of the spectrum is likely to occur, we have found that fits 
to data weighted by a single factor of k give more reliable 
results overall. In general, if we are not focusing on the 
0(4) part of the spectrum, we fit our data to the FT of 
kx(k). Fittiag to the FT reduces the effect of overlap- 
ping peaks £3 Since XAFS can be thought of as a sum of 
the contribution to the absorption of the individual atom 
pairs, we calculate F(k), A(fc), S c (k) and Si(k) given an 
approximate atom cluster for each scattering path us- 
ing FEFF6. This information for each path is called a 
"standard." The data are fit to a sum of theoretically 
calculated standards which include all single scattering 
paths up to the Cu(l)-Cu(2). The only multiple scat- 
tering paths that are included involve the Cu(l)-Cu(2) 
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path and the intermediate 0(4). The Cu(2)-planar oxy- 
gen path between the planes was too low in amplitude 
to obtain reliable fit results, so its parameters were fixed. 
Other multiple scattering paths were found to be negli- 
gible, or were dominanted by the single scattering paths. 

Several of the parameters were constrained in order to 
decrease the number of fit parameters. The Cu(2)-Y and 
the Cu(2)-Cu(2) distances were constrained such that the 
yttriums rest in the center of the four Cu(2) neighbors 
given the diffraction measurements of the a and b lat- 
tice parameters of typical samples. This constraint was 
necessary to fit the Cu(2)-Cu(2) peak (3.38 A), which is 
small compared to the nearby Cu(2)-Y (3.21 A), Cu(2)- 
Ba (3.37 A), and Cu(l)-Ba peak (3.48 A). Also, if the Cu- 
Ba peaks were allowed to vary over a broad range, they 
would often fall into a false minimum in the fitting pa- 
rameter such that they were ~ 0_04 A too long. This be- 
havior was noted by Stern et aln However, we were able 
to determine by fits to a Cul standard (a Cu-I standard 
should be very similar to a Cu-Ba one) that this error 
was due primarily to the overlap of several neighboring 
peaks in YBCO, and not to unusual errors in the FEFF6- 
generated backscattering amplitude. We therefore only 
allowed the Cu-Ba peaks to vary over a narrow range 
(±0.02 A from the diffraction result). All width param- 
eters were allowed to vary freely, except o"c«(2)-o(2,3) 
(across the planes), which was held fixed. Eq's for like- 
atoms were held equal, but were allowed to vary between 
different types of atoms to help correct for minor errors 
in the FEFF6 phase calculations. 

The aim of these fits is to identify parameters that 
change with temperature. Since measurements of Sq and 
a are correlated, we performed a fit where these parame- 
ters were allowed to vary freely, and then calculated the 
average for Sq for each peak. A similar correlation occurs 
between measurements of Ri and Eq, so average Eq's for 
the individual peaks were also determined. The values 
of Sq and Eq were then fixed for each path and the fits 
were performed again. By holding Sq and Eq fixed for all 
temperatures, we ascribe all changes in the peaks' am- 
plitude and frequency with temperature to changes in cji 
and R. 

Although Sq should be the same for all paths in prin- 
ciple, we found that if a single value of Sq was used 
high quality fits could not be obtained for the Cu(l)- 
Cu(2) peak (which includes multiple scattering off the 
0(4) atom). The fit was improved by allowing this one 
peak to have an Sq=0.75, while all other peaks use an 
Sg=0.90. 

The maximum number of parameters given our fit 
range using the method of SternB is 29. The total num- 
ber of parameters in these fits was 17, well below the 
maximum. 

Errors in the parameter measurements are difficult to 
estimate reliably in XAFS fits. We estimate the errors 
on all the parameters from the covariance matrix gen- 
erated by the fit. The variance of the data used in the 
covariance matrix is obtained by assuming that the resid- 



ual difference between the data and the fit is normally 
distributed. This procedure only accounts for random 
errors and thus should describe relative errors between 
data at nearby temperatures. Systematic errors due to 
problems in the FEFF6 calculation, self-absorption cor- 
rections, estimates of /io, etc., will cause overall shifts in 
the best fit. Absolute errors on these measurements are 
generally < ±0.02 A in R and < ±15% io-Sq and a, and 
are generally better for nearest neighbor.EiJ 

B. Fitting results 

1. 0(4)site distribution in the single crystal 

As suggested by the presence of the beats, a two 0(4) 
site distribution was necessary to fit the single crystal 
data in Fig. ||for the four lowest temperatures. The fits 
shown in Fig. || all assume a two-site distribution for 
the 0(4) except for the T=100 K data. For the T=100 
K data, a single site for the 0(4) produces a reasonable 
fit; the slight dip in the amplitude at ~ 11 A^ 1 can be 
modeled by an interference between a single Cu(l)-0(4) 
path and a single Cu(2)-0(4) path. Note that the data 
in this region is still visibly different from the film data. 
The features in the T=50, 78, 83, and 86 K data between 
10 and 13 A" 1 are modeled by two Cu(l)-0(4) paths and 
two Cu(2)-0(4) paths, with each path length split by 0.12 
± 0.01 A and 0.06 ± 0.06 A, respectively. We did not 
observe any obvious trend with temperature other than 
for the T=100 K data. 

The two-site fits to the single crystal data may be de- 
ceptively good. In order to obtain these fits, unphysically 
small (< 0.01 A) values for the broadening factors were 
needed. In addition, the sum of the amplitudes for these 
peaks (which should add up to the number of neighbors 
times Sq) is about 50% too high. These exceptionally 
narrow peaks indicate that besides the split peak distri- 
bution, further anharmonicity must exist. A common 
way to model anharmonicity in XAFS data analysis is to 
expand Eq. ^ about Ri for a given atom-pair in pow- 
ers of < r" >, otherwise known as a "cumulant expan- 
sion" (< r™ > is the n th moment, M4 is the part of the 
4th moment that is_different from a gaussian, namely 
M A =< r" > -3cr 4 )Jd The 2nd cumulant is the width in 
the harmonic approximation, or the Debye- Waller factor. 
The 3 rd cumulant affects the phase shift. The 4 th cumu- 
lant multiplies the overall amplitude by a factor es k Mi . 
If this term is positive, it increase the amplitude of the 
XAFS at higher wave vectors, acting like an imaginary 
Debye- Waller factor, or like a distribution with a cusp 
which puts more weight at the center of the distribution 
at the expense of the wings.Ej Excellent fits to the single 
crystal data can be obtained by allowing the four Cu- 
0(4) peaks to have Af 4 = 2.7 x 10 -5 A 4 . These fits have 
much more reasonable Debye- Waller factors (^0.04 A) 
and Sq is no longer 50% too large. With such a large 
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value of M4, the cumulant expansion is no longer accu- 
rate above a k of ~ 14 A. 

These fits serve to illustrate that we can model the 
0(4) site distribution with reasonable values for the fit 
parameters, but that the fits are not unique. Because Sq, 
a and < r A > are so correlated for these slightly separated 
Cu-0(4) peaks, a study of the temperature dependence 
of any of these parameters would produce questionable 
results. The only firm result is that at least three, and 
possibly four, atom-pair distances are necessary to de- 
scribe the beat structure in the /c-space transforms for 
the single crystal. The actual 0(4)-site distribution g(r) 
may be more complicated. 

Table | reports fit results to the further neighbors for 
the T=50 K sample. 



2. Fits to the film on MgO 

The Cu(l)-0(4) and Cu(2)-0(4) peaks are well de- 
scribed by a single 0(4) site in the film data, as suggested 
in Sec. |V] (Fig. M). Allowing for a two site distribution 
with a A r > 0.09 A did not produce reasonable fits. Fits 
to the film data of a similar quality as to the single crys- 
tal data are obtained without assuming any anharmonic 
behavior in the atom-pair distances out to ~ 4.1 A. Fig. 
H shows the fit to the film data over this range both in 
fc-space and in r-space, and Table | reports the fit to all 
the parameters for both the single crystal and the film 
on MgO at T=50 K. No significant deviation between 
the XAFS atom-pair distance measurements and stan- 
dard diffraction measurements were observed, except for 
the unusually short Cu(2)-0(4) and Cu(l)-Cu(2) pairs, 
which are about 0.02 A shorter than in the single crystal 
diffraction measurements. This measurement is consis- 
tent with the shorter c-axis of this material (11.64 A at 
room temperature) compared to the diffraction sample 
(11.68 A at room temperature). 

Plots of a vs. T are shown in Fig. ^. All the atom- 
pairs show an increase of a with temperature, except 
the Cu(l)-0(4) and the Cu(2)-Cu(2) (not shown) atom 
pairs. The measurements of 0"c?«(2)— C«(2) are fairly large 
(~ 0.11 A) with large estimated errors and are there- 
fore unreliable. Only the Cu(l)-0(4) and the Cu(2)-0(4) 
pairs show any anomalies in a near T c . a for both pairs 
increases slightly just above T c and then returns just be- 
low T c . The jump is larger for the Cu(2)-0(4) width 
(~0.008 A) than for the Cu(l)-0(4) width (~ 0.004 A). 
This behavioc.was noted for the Gu(l)-0(4) bond both 
by Stern et au and in Kimura et a£j, although the Stern 
result assumed a split 0(4) site. None of the previous 
XAFS studies report a change in the Cu(2)-0(4) width. 

All other atom-pair parameters show smooth behavior 
with temperature and are of reasonable values. 



C. Correlations of the further neighbors 

We have analyzed the correlations between positions of 
the Cu atoms and their near-neighbors by calculating </> 
for each atom pair. Diffraction measurements of broad- 
ening factors for each site are required, as indicated in 
Eq. We have chosen temperature-dependent neutron 
diffraction data given by Sharma et a£ZI because they in- 
clude a relatively dense grid in temperature between data 
points (AT =10 K), measurements over a similar range of 
temperatures (10-300 K) to this study, anisotropic ther- 
mal factors for the 0(4) site, and a high sensitivity to the 
oxygen atoms. This diffraction study saw no anomalous 
features with temperature in the thermal factors for any 
lattice,site and is consistent with a similar study by Kwei 
et alEB Since Sharma et aPs measurements of the thermal 
factors were not always at the same temperatures as the 
measurements in this report, we fit each Debye- Waller 
factor from diffraction vs. temperature with a polyno- 
mial. The use of these thermal factors may introduce 
systematic errors into our measurements of <j> not only be- 
cause the samples are prepared with different methods in 
different laboratories, but also because the diffraction re- 
sults were obtained from a powder of YBCO and our data 
is for a film. Unfortunately, we must tolerate such er- 
rors because no comprehensive, temperature-dependent 
study of a thin film of YBCO exists at this time. Such 
errors should only contribute to an overall shift of the 4> 
parameter and should not greatly affect the temperature 
dependence unless major features in the (as yet, unmea- 
sured) diffraction thermal factors occur in measurements 
on films and not on single crystals. 

The correlation coefficients vs. temperature for each 
atom pair are displayed in Fig. [l^. Errors are propa- 
gated from the quoted errors in Ref. ^ and the errors 
in Fig. ^ and in no way attempt to reflect any system- 
atic errors introduced by FEFF6, differences in the sam- 
ples, etc. The correlations for the Cu(2)-Cu(2) bond are 
not displayed because our measurements of 0"cu(2)— Cu(2) 
are unreliably large with large estimated errors. All the 
other atom-pairs are resolved in the fits well enough to 
give reasonable estimates of <f). 

VI. DISCUSSION 

A. Split 0(4) site 

The main result reported above is that while two 0(4) 
sites are necessary to describe the XAFS data for a sin- 
gle crystal very well, a high-quality fit to XAFS data for 
the thin film with a single 0(4) site is also possible. Al- 
though this result only puts an upper limit on the size of 
any site splitting at <0.09 A, the measurements of the 
Debye- Waller factors are small enough to indicate that 
any possible splitting is probably even smaller. (Any un- 
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resolved splitting should add to the Debye- Waller factors 
in quadrature.) 

Even though the existence of the splitting in some sam- 
ples is now well established, the nature of the splitting 
is still quite unclear. The main problem is that the split 
sites could be due to motions of individual atoms between 
distinct sites, or due to two separate harmonic potentials, 
displaced fcam each other by, perhaps, a local distortion. 
One modelli3 for static splitting relies on oxygen vacancy 
ordering on the 0(1) site, which may shift the local po- 
sition of the 0(4) atom. At least one configuration of 
oxygen vacancies can produce a split 0(4) position when 
half of the 0(1) sites are occupied: the Ortho-II phase of 
YBa 2 Cu 3 06.5, which occurs when every other Cu-O(l) 
chain is completely vacant of oxygen. In this case, only 
half the 0(4)'s have near 0(1) neighbors. This config- 
uration splits the 0(4) site into two sites separated by 
±0. 05.-4 from_the stoichiometric (no 0(1) vacancies) 
site.E£fE3 RohlerEj considers a similar situation when sin- 
gle 0(1) vacancies cause shifts in the 0(4) position both 
along the same chain and along neighboring chains. In 
this model a single 0(1) vacancy affects 12 0(4) atoms. 
In the case of optimally doped YBCO with S = 6.93, this 
would mean 84% of the 0(4) oxygens could be affected. 
This model also helps explain why the effect seems to 
be more easily fit as a split 0(4) distribution in lower 
T c sampleSp/i.e. more 0(1) vacancies) as measured by 
Stern et aln Any changes observed in the splitting near 
T c would be considered a reordering of oxygen vacancies 
in this model. In order to explain the thin film results, 
one could invoke the added stress on the material intro- 
duced by the substrate. This stress could perhaps be 
manifest as a freezing of the vacancies into a random or- 
ganization that suppresses the amount of splitting. This 
hypothesis is not likely, however, because oxygen diffu- 
sion measurements aife,pot drastically different in a film 
compared to crystals.LJO Another possibility is that the 
films are overdoped and hence there are very few 0(1) 
vacancies. A direct test of this argument would be to 
measure films that range from overdoped to underdoped. 

A dynamic model that would lead to either a broad- 
ening or a split 0(4) position is the hopping small po- 
laron model, with the polaron hopping on and off the 
Cu-0(4) bonds as discussed by Mustre de Leon et alH 
Ranningero has calculated the XAFS signatures for the 
presence of polarons and shown that as the hopping speed 
decreases, the pair distribution function first broadens 
and then splits; less than a factor of two in hopping speed 
is required to change from a broadened peak to a split 
peak. This is expected from the following simple argu- 
ment. When the electron hops at a much faster rate than 
the lightest optical phonons, the lattice cannot respond, 
and there is little displacement of the individual atoms. 
If the hopping is slower, the lattice has time to respond, 
and the Cu and 0(4) atoms move towards each other or 
apart as the polaron hops on and off the ligand. If the 
hopping rate falls below the optical frequency, a split Cu- 
0(4) distribution emerges, corresponding to the absence 



or presence of the polaron. This model can easily explain 
both the single crystal and the thin film results with a 
small change in the polaron hopping rate. The fact that 
the Cu(2)-0(4) peak is less ordered than the Cu(l)-0(4) 
peak and has a correlation parameter ^0.5 is consistent 
with such a model. 

Although the oxygen vacancy argument has possible 
merit, there is a body of literature that links dynamics 
of the 0(4) sitaln-anDmalies in vibrationaLsjpejctrpacopies 
both in YBC0t3BE3 and other material. ESfelBEfl Since 
XAFS is not sensitive to static verses dynamic changes 
in atomic positions, there is little chance the XAFS tech- 
nique can resolve this issue. However, XAFS can deter- 
mine whether a pair-distribution function is harmonic or 
not. In our fits to the split distribution in the single crys- 
tal, we had to allow for a large fourth cumulant for both 
0(4) sites. This term can either have the effect of flatten- 
ing out or sharpening the peak of an otherwise harmonic 
distribution, depending on the sign. These data required 
a positive 4th-cumulant, and were thus sharper than a 
harmonic distribution. No third cumulant was necessary 
to fit the data. The necessity of this extra parameter in- 
dicates that the pair-distribution function is not simply 
the sum of two harmonic distributions. The need for this 
quartic term is consistent with Raman measurements of 
the 0(4) mode at 505 pip" 1 which also requires a sig- 
nificant 4th=prder term.£3 This anharmonicity has been 
speculated!!^ to be relatecL|to the anharmonic potential 
implied by the split 0(4).Bc3 The anharmonicity in the 
0(4) 505 cm -1 mode could now be directly linked to the 
split 0(4) distribution if accurate measurements of the 
Raman anharmonicity can be shown to be significantly 
different between a film and a single crystal. This result 
would mean that the anharmonicity (and any polaron for- 
mation related to this anharmonicity) is not required for 
high-T c superconductivity. On the other hand, if the an- 
harmonicity in the 505 cm" 1 peak is unchanged between 
the film and the crystal, then the split 0(4) is not related 
to it. In this scenario, polaron formation may still exist 
and be important for high-T c , but the split 0(4) would 
not be the structural manifestation of it. 

When trying to explain why splitting of the 0(4) site 
may or may not be present, one should also consider the 
T c dependence. Stern et al reports that the changes in 
the splitting near T c are more pronounced for the lower 
T c samples. However, the two-site fits to the 0(4) for 
the higher T c samples in their work are not of the same 
high quality as for the lower T c samples. In other words, 
the pair-distribution function for these samples appears 
to be more complicated than a simple two-site distribu- 
tion. Therefore, the Stern fits cannot really rule out any 
significant temperature dependence. Certainly our own 
results on the single crystal indicate that the tempera- 
ture dependence may not be consistent from sample to 
sample; all the data at T < T c cannot be fit satisfactorily 
with a single 0(4) site, yet the T=100 K data is well de- 
scribed with a single 0(4) site. This result is in contrast 
to the work by Mustre de Leon et al, where they report 
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the splitting to be constant at all temperatures except 
near T c where the splitting shrinks from 0.13 A to < 
0.11 A, and to the work by Stern et al, where for some 
samples, they report a constant splitting at all tempera- 
tures. 



B. Changes in broadening parameters 

In our measurements of the broadening parameters of 
the film on MgO we observe some significant temperature 
dependencies. Almost all the measured bonds show the 
Debye- Waller factor increasing with temperature, as one 
expects if the Debye temperature is not too large (Fig. ||). 
The exception is the Cu(l)-0(4) bond which maintains 
a broadening of approximately 0.035-0.40 A from T=20 
to 200 K. The only significant behavior occurring near 
T c involves the 0(4) atom: both the Cu(l)-0(4) and the 
Cu(2)-0(4) pairs show an increase in their broadening 
parameter, with ccu(i)-0(4) jumping from 0.035 to 0.040 
between 80 K and 100 K and back then to 0.035 A, and 
CT Cu(2)-0(4) jumping from 0.060 A to 0.068 A between 80 
K and 96 K and back to 0.060 A. In the region between 
T=80 K and 100 K, both (t Cu (i)~o(a) and <7 Cu{2 )-o(A) 
varies from temperature to temperature more than the 
estimated error. This fluctuation may be real or the er- 
ror estimates may be too small for some reason, perhaps 
from inaccurate temperature measurements. However, 
the Cu(2)-0(4) data vary more smoothly. Stern et al 
also measured a fluctuation in the Cu(l)-0(4) broaden- 
ing near T c , however their measurement assumed a two 
site distribution that had collapsed to one site in the 
same temper.ature range as the broadening fluctuation. 
Kimura et a££l measured a fluctuation in <7cwi)-o(4) m 
a similar temperature range using a single 0(4) site fit, 
but reported no anomalies in <Jcu{2)-o(i) near T c . No 
proven explanation of the physics behind these fluctua- 
tions exists, however, they can be interpreted in terms of 
changes in the correlation between the Cu and O atomic 
positions (see next section). 



C. Correlations between atomic positions 

Like the 0(4)-site position, changes in <j> with temper- 
ature can be attributed to changes in the static disorder 
of the individual sites, or to changes in the dynamics 
of the lattice displacements. As the temperature of the 
sample is decreased, static disorder could manifest itself 
as an ordering of oxygen vacancies, which would then 
generate a distinct set of positions for a given atom. If 
the site is split into two or three distinct positions, XAFS 
would measure an unusually broad peak for that pair and 
thus look like the atom-pair displacements are negatively 
correlated. It is therefore important to consider the ab- 
solute value of the broadening measurements from XAFS 
as compared to diffraction to try to help determine if the 



model used to fit the data is consistent, i.e. single atomic 
sites broadened harmonically around some average pair 
distance. Also, trends in the correlation with temper- 
ature should give some insight into which correlations 
within the unit cell are dependent on static or thermal 
disorder. 

The most important result from the measurements of 
<fi in Fig. [l^ is that, within the estimated errors, they are 
all in the range from to 1, indicating that the measure- 
ments of Sq, the calculated F(k) and the model used to 
fit the data are all fairly accurate. More direct tests of 
the reliability of the absolute measurements of <f> are pro- 
vided by the nearest and the furthest bonds measured, 
that is, the Cu(l)-0(4) pairs and the Cu(l)-Cu(2) mul- 
tiple scattering pair. Since the Cu(l)-0(4) pair is the 
nearest neighbor pair in YBCO, we expect at low tem- 
peratures that 4> should be very near unity, and it is: 
4>Cu{X)— 0(4) demonstrates weak temperature dependence 
with a mean value of about 0.87, in approximate agree- 
ment with correlation measurements of the Hg-0(2) pair 
in Hg-120lQ The fluctuations near T c are well within the 
increased errors after propagating the diffraction errors 
into <t>. 

The Cu(l)-Cu(2) atom pair is the furthest neighbor fit 
and includes multiple scattering off the intervening 0(4) 
atom. It is therefore a good test of the FEFF multiple 
scattering calculations for F(fc)c3 as well as the absolute 
reliability of </>. We expect this pair to be the least cor- 
related, and it is; 4>cu(i)-Cu(2) shows no obvious change 
with temperature (although both XAFS and diffraction 
broadening factors show a significant temperature depen- 
dence) and is consistent with 4>Cu(i)-Cu(2) — 0. The ab- 
solute accuracy of this pair is worse than the other pairs 
measured because the FEFF calculation includes multi- 
ple scattering off the 0(4). Consequently, this measure- 



ment is consistent with . 



>Cu(l)-Cu(2) 



S ±0.2. The Cu-Ba 



pairs show a strong temperature dependence in <fr, in each 
case dropping from about <fi=0.7 to 0.5 between 20 and 
200 K. The Cu(2)-Y pair is partially correlated (~ 0.45 
±0.1) and shows little temperature dependence between 
20-200 K. The lower correlation coefficient for the cold 
temperature Cu(2)-Y pair compared to the Cu-Ba pairs 
is probably because both atoms in the pair are bound to 
the planar oxygens, with bond angles between the Cu(2)- 
0(2,3) and Y-0(2,3) bonds ^90°; the Ba atoms are more 
constrained by the axial oxygens. 

The most interesting correlation measurements are for 
the Cu(2)-0(4) pair. 0c«(2)-O(4) starts out very low 
(~0.4) but starts to increase to a maximum of ~0.6 just 
below T c . In the vicinity of T c , (f> then drops to about 
0.45 and then increases again to its maximum value of 
0.6. As the temperature is increased above 100 K, 4> de- 
creases steadily from 0.6 to 0.4. This decr ease i n <f> (and 
hence the increase in a mentioned in Sec. VI B) may be 
the result of a negatively-correlated mode being excited 
near T c relative to the correlated mode that is dominat- 
ing away from T c . Such a mode mav-be consistent with 
a polaronic-hopping transport modea 







VII. CONCLUSIONS 

We have contrasted XAFS measurements on a thin film 
of YBCO on MgO with measurements on a single crystal 
with 2% Ni. We have also compared the results with a 
film on LaA103, which gives essentially the same results, 
and with previous measurements on oriented powders. 
The single crystal data exhibit an anharmonic site dis- 
tribution for the 0(4) atom which can be described as 
a split position below T c and a single position above T c . 
(The temperature range measured for the single crystal 
is limited.) This result-is consistent with previous results 
on oriented powders .Em 

The main result of this paper is that the thin film data 
show very different behavior. The 0(4) peak can be well 
described by a single-site, harmonic distribution. This 
distribution does not show any strongly anharmonic be- 
havior with temperature near T c , or at any other tem- 
perature between 20-200 K. 

These data refute the argument that the XAFS signal 
showing the split 0(41 site distribution is too small to 
be measured reliablyp however, high S/N is required. 
The Fourier-filtered film data are very reproducible over 
the temperature range measured, and since only subtle 
changes occur in the 0(4) site distribution for the film 
near T c , the changes in the data with temperature are 
gradual and predominantly thermally driven. Because of 
this reproducibility, the rapidly changing XAFS from the 
single crystal must be taken as a legitimate signal, and 
not as noise. 

Although we don't see any evidence of anharmonicity 
in the 0(4) site in the film, fits show that the broadening 
factor for the Cu(2)-0(4) bond exhibits a local maxi- 
mum above T c followed by a drop below T c to essentially 
its T=20 K value. The Cu(l)-0(4) pair demonstrates 
similar behavior, but the effect is smaller and is not as 
systematic with temperature. Since there is no .such be- 
havior in the published diffraction literatureEZrEa we have 
interpreted this as a change in the degree of correlation 
between the Cu(2) and the 0(4) atoms, and possibly be- 
tween the Cu(l) and the 0(4) atoms. This result is in 
contrast to measurements by Kimura et a££l which show 
changes in the broadening factor for the Cu(l)-0(4) bond 
near T c , but not in the Cu(2)-0(4) bond. 

Finally, a better understanding of the local structure 
is obtained by considering the local correlations of the 
near-neighbors. As mentioned above, the Cu(l)-0(4) is 
measured to be a very tight bond, yet the Cu(2)-0(4) 
is "looser" and exhibits a decrease in the correlation be- 
tween the atoms near T c , indicating that the chains act as 
a unit that is somewhat independent of the planes. This 
lack of correlation is also indicated by the measurement 
of the correlation coefficient 4>cu(i)-Cu(2) near zero (< 
0.2) that also does not change with temperature. How- 
ever, the atomic displacements in the metal-oxide layer 
do become less correlated at higher temperatures. The 
atomic positions in the Cu02 planes apparently remain 



correlated in the temperature range measure based on 
the measurements on the Cu(2)-Y atom pair. 
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TABLE I. Comparison between fit results to the single crystal of YBCO:Ni 2% (xtal) and the filrn-ef YBCO on MgO 
at T=50 K. n.b. stands for "number of bonds per unit cell." Neutron diffraction results of Sharm et a£3 are also given for 
comparison, and measurements of <f> are reported for the film data. Measurements of 4> for the single crystal are qualitatively the 
same. Errors marked with an "x" are unreliable because the pair-distribution function deviates significantly from a gaussian. 
Imaginary Debye- Waller factors often accompany these distributions; such parameters caus e the XAFS oscillations to increase 
with k (unphysical). Some parameters were held fixed or constrained in the fits. See Sec. VB for further discussion. Direct 



comparisons of bond lengths should consider that at room temperature, the diffraction sample has a c-axis lattice parameter 
of 11.68 A, while the film's is 11.64 A. The lattice parameters of the single crystal are not known. 
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FIG. 1. The crystal structure for YBa— 2Cu307-i . The 0(2) and 0(3) oxygens are often referred to as the "in-plane" 
oxygens. Much of this paper is concerned with the 0(4), otherwise known as the "axial" oxygen. 
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FIG. 2. Full e-space data for a YBCO film on MgO at T=20 K with the x-ray polarization parallel to the c-axis of the film. 
These data have been dead-time corrected and had some minor glitches removed, as described in the text. Data from the film 
on LaAlC>3 and the single crystal are of comparable quality. 
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FIG. 4. The Fourier transform (FT) of k 3 X (k) (= X 3 M) vs. r for (a) the same sample as in Fig. ^, and (b) the single 
crystal of YBCO: Ni 2%. Data are transformed with a Gaussian window between 3-15.5 A -1 and broadened by an additional 
0.3 A -1 . The oscillating curve is the Re(x 3 ), and the envelope of this curve is the amplitude, [Im(x 3 (r)) 2 + Re(x 3 (r)) 2 ]z . 
The peak at 1.55 A corresponds to the Cu(l)-0(4) atom pair, while the shoulder of this peak at 2.0 A corresponds to the 
Cu(2)-0(4) atom pair. The multi-peak between 2.4 and 3.4 A is a combination of the Cu(2)-Y, Cu(2)-Cu(2), Cu-Ba and 
Cu(2)-0(2,3) atom pairs. Lastly, the peak at 3.8 A is the Cu(l)-0(4)-Cu(2) multiple scattering peak. The peak positions are 
shifted from the actual pair-distances because of the combined effect of 5 C , Si and Fi(k) in Eq. El 



1G 




FIG. 5. 0(4) contribution to k 3 \{k) vs. k at several temperatures for the single crystal of YBCO: Ni 2%. The solid lines are 
the data and the dotted lines are the fits. Each fit assumes a split 0(4) site except the T=100 K fit. The 0(4) contribution 
was obtained in an identical fashion to the data in Fig. H. 
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FIG. 6. 0(4) contribution to k 3 \{k) vs. k at several temperatures for the YBCO film on MgO. The 0(4) contribution was 
obtained by back-transforming r-space data (Fig. ^) between 1.3-2.2 A. Contributions from high shells (say, from the Cu(2)-Y 
peak) have been shown to be negligible.!!! Only half the data is displayed for clarity; an additional scan was taken in between 
each of the temperatures displayed here. 
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FIG. 7. Fit to the 0(4) contribution to k 3 X (k) vs. k for the YBCO film on MgO at T=20 K. Fit assumes a single 0(4) site. 



19 




20 



0.045 - 



0.040 • 

S 



0.035 - 



Cu(l)-0(4) 



s 



<) 

-- <l> 



0.030 — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i— i — i — i — p 

50 100 150 200 

r(K) 



0.08- 



0.07- 



0.06- 



0.05 



Cu(2)-0(4) 0o 



(I (I 



i — i — i — i — | — i — i — n — | — i — i — i — i — | — n — i — i — | — r 
50 100 150 200 

T(K) 



0.055 - 
0.050 
C 0.045 
0.040 
0.035 H 



Cu(2)-Y 



o o l- 1 - 



l 1 — I — i — I — I — I — I — i — I — I — I — I — i — I — I — I — I — i — r 

50 100 150 200 

r(K) 




t — i — i — i — i — i — i — n — i — i — i — i — i — i — n — i — i — | — r - 

50 100 150 200 

r(K) 



0.07- 



0.06- 



v 0.05 - 



0.04- 



0.03- 




i — i — i — i — | — i — i — n — | — i — i — i — i — | — n — i — i — | — i 
50 100 150 200 

T(K) 



0.08- 



0.07- 



0.06 




i — i — i — i — | — i — i — i i — | — i — i — i — i — | — i i — i — i — | — i 

50 100 150 200 

T(K) 



FIG. 9. a vs. 

Cu(2)-Cu(2) (at r 



temperature from the film of YBCO on MgO for all single scattering paths up to 4.15 A, except the 
-3.38 A) and the Cu(2)-0(2,3) (at -3.66 A]_ pairs. Errors are estimated from the covariance matrix, and 



in no way reflect any systematic errors, as discussed in Sec. 
indicate the reproducability of the data. 



VA 



These errors should thus be taken as relative errors which 
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FIG. 10. The correlation coefficient 6 vs. temperature from the film of YBCO on MgO for the same atom pairs as in Fig 
Errors are propagated through Eq. g using the errors for the diffraction data reported in Ref. |37]. 



22 



